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Abstract 

The inverse scattering theory for the sine-Gordon equation discretized in space and 
both in space and time is considered. 

1 Introduction 

In the framework of integrable nonlinear evolution equations, Hirota recovered almost 25 
years ago the Bianchi superposition formula for the solutions of the sine-Gordon equation, 
showing that it can be considered as an integrable doubly discrete (both in space and time) 
version of the sine-Gordon equation He found also its Lax pair, Backlund transformations 
and 7V-soliton solutions. In the following, this study was generalized to include the discrete 
(in space or time) case soliton solutions were obtained by using the dressing method 
, the related inverse spectral transform was studied and the nature of the numerical 
instabilities of the solutions was discussed ||, |^ . 

Here we are interested in the spectral transform for the discrete and doubly discrete 
sine-Gordon equation with solutions decaying to (modvr) at space infinity . 

The aim of this paper is twofold. 

From one side, we re-examine the scattering theory of the spectral problem introduced 
by Orfanidis in ||^ which depends on a discretized space variable n. According to the usual 
scheme, we introduce two summation equations which define a couple of Jost solutions /in 
and Vn characterized by their asymptotic behaviour, respectively, at n — > —oo and n ^ +oo. 
These summation equations are, as usual, used to define the spectral data. However, we 
show that, in contrast with the continuous case, in order to study the analytical properties 
of the Jost solution it is necessary to introduce an additional summation equation that 
only in the continuous limit reduces to the previous one. Then, we find the conditions on 
the potential that ensure the existence and analyticity properties of the Jost solutions, we 
formulate the inverse problem as a Riemann-Hilbert boundary value problem on the real 
axis and we derive the time evolution of the spectral data, paying special attention to the 
nature of the singularity of the time evolution of the solutions a.t t — 0. 

From the other side, we show that both in the semi-discrete and doubly-discrete case, the 
equations of the matrix Orfanidis spectral problem can be decoupled and reduced to a scalar 
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problem, giving as principal spectral problem the "exact discretization" of the Schrodinger 
operator introduced by Shabat Q 

</'«+2 = 5«'/5„+i - (1 + A)(p„ (1.1) 

and obtained iterating the Darboux transformations. Then, following a procedure analogous 
to that used in |^ for getting a 2+1 generalization of the sine-Gordon equations, we do not 
consider a Lax pair but a Lax triplet, i.e. we derive the discrete sine-Gordon equations as a 



compatibility condition of (|l . l| ) with a pair of auxiliary spectral problems 



The spectral theory for the operator (LI) was exhaustively studied in [g| for a potential 
(7„, not necessarily real, satisfying 

-|-oo 

(l + |n|)|g„_i-2| <«). (1.2) 

n— — oo 

and applied to a discrete version of the KdV. However, one must be advised that the solutions 
of the discrete sine-Gordon equations are related to the potential gn via an equation that 
can be considered a discretized version of a Riccati equation and can be solved in terms of a 
continous fraction. Moreover, the potential g„ is complex and the problem of characterizing 
the spectral data is left open. 

By using this theory we introduce the spectral data, we find their time evolution and, 
therefore, according to the usual inverse scattering scheme, the Cauchy initial value problem 
is linearized. In particular, we show that the time evolution is discontinuous at the initial 
time t = 0. 

We are also able to find the couple of doubly discrete auxiliary spectral operators that in 
triplet with the doubly discrete Schrodinger operator gives the Hirota-Bianchi fourth order 
doubly discrete sine-Gordon equation. 

The whole theory can be trivially extended to the sinh-Gordon case. Then, the potential 
Qn is real and the characterization equations for the spectral data are easily obtained in 
analogy with the continuous case. 

2 Orfanidis Lax pair 

Let us consider the Lax pair proposed by Orfanidis |^ 

Xn,m-|-1 — ^n,m\n,7n (^■^) 



with 



g-»(e„+i,,„-e„,„)/2 

ik e'(^"+i--''-")/2 



(2.3) 



/ 1 Tlg-i(0.>.™ + l+0„.,„)/2 \ 

\ ik J 

with 0„,m real, k the spectral parameter and 7 a real constant. The compatibility condition 
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where E shifts n and F shifts m gives by inserting in it (2.1) and (2.2) Mn,m+iNn,m ~ 
Nn+i,mMn_m- From it we have the Bianchi-Hirota equation 

/ ^n+l.m+l ^n+l.m ^n,m+l ^" ^n,m \ . / ^'n + l.rn+l ^ ^n+l.m ~t~ ^n,m 

sm ^ — — 7 sm 



(2.5) 

If we introduce 

e,^j = e,{t + Tj), t = TOT, (2.6) 

change the constant 7 as foUows 

7 ^ T7 (2.7) 

and take the Umit r ^ and ni 00 &i t fixed we get the semi-discrete sine-Gordon 
equation in the Ught cone coordinates 

dtOn+i - dte^ = 47 sin i (0„+i + 0„) . (2.8) 
Notice that this equation can be obtained as compatibility condition of the Lax pair 

Xn+l = MnXn (2.9) 

dtXn = N^Xn (2.10) 

with 

^"=( .fc e^(.^.^-M/0 (2-11) 

Notice also th at in both cases, doubly discrete and space discrete cases, the principal 
spectral problem ( ^.l| ) is the same. 

3 Direct scattering problem 

3.1 Jost solutions and summation equations 



Let us rewrite the spectral problem ( 2.1 



Xn+l - + ik(Ti)Xn = QnXn (3.1) 

where the Pauli cr matrices are defines as usual, i.e. 

-(? «")• "-(J -0 

and 

Q„ = e-3'^='(»"+i-»") - 1 (3.3) 
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and let us consider it for k = kj^c- 

If Gn is a Green function, i.e. if G„ satisfies 



Gn+l - {I + tkc7i)Gn ^ dnfll (3.4) 



a matrix solution of (3.1) is given by the solution of the following summation equation 

+00 

Xn=Wn+ ^ Gn-jQjX] (3-5) 
j = -oo 

where Wn is a solution of the homogeneous equation 

Wn+i - (1 + ikai)wn = 0. (3.6) 
We represent G„ and 6nfi as Fourier integrals 



Gn = ^i p^-'G{p)dp, (3.7) 
27ri J\p\=R 

<o ^7^<f P"''dp. (3.8) 

2tti J\p\=R 



Inserting into (p^) we get 



and therefore 



{{p - 1)1 - ikai)G{p) ^ 1 (3.9) 



G{P)= ^^_^^,^^, ((p-l)l + »fcai). (3.10) 
Notice that G{p) has a pole dX p = \ ±ik and notice also that at large p 

l+zfccTi y^ (l + zfccTl)^ ^^^^^ 



and therefore at large p 



(3.12) 



We consider k — kjic and for different choices of R, i.e. for R > 1 and R < 1, we get 
different Green functions, G and 7?, and, correspondingly, different summation equations 
defining different Jost solutions, (f) and ip. 

For R > 1 + fc^ the integrand in (3.7) inside the disk |p| < i? contains a pole at p = 
for n < and the poles p — 1 ±ik and it is analytic outside the disk. The integral can be 
computed evaluating the residuum at p = oo by using ( 3.12 ) and we get 

Gn = Q{n - I) {1 + ikaif^ (3.13) 

where &{n) is the discrete version of the Heaviside function 
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Therefore 

n-l 

<j)n = Wn+ ^ {l + ikai)"~^~^ Qj4ij. (3.15) 

j=-oo 

If < 1 we have to subtract to the prevfous integral the residua at the poles p = 1 ± ik. 
We have 

Hr, = e(n - 1) (1 + ikaif-' - ^(1 + ai)(l + ik)^-^ - ^(1 - ai)(l - ifc)"-^ (3.16) 



Since 



we have 



(1 + ikaif = i(l + ai)(l + ik)" + i(l - ai)(l - ifc)" (3.17) 



Hn = -e(-n) (1 + ifccTi)""^ (3.18) 
and the summation equation becomes 

^n = Wn-^il + ikuiy-'-^ Qj^j. (3.19) 

It is convenient to choose 

u;„ = (l-iCT2)(l + ifc(T3)" (3.20) 

so that (j)n and tpn satisfy the following boundary conditions 

^„ ~ (1 — i(T2) (1 + jfccTs)" n — > —00 (3-21) 

ipn^i^- i(T2) (1 + ifccrs)" n ->■ +00. (3.22) 

In order to study the analytical properties of the Jost solutions with respect to the spectral 
parameter k, it is necessary to introduce the modified matrix Jost solutions 

Hn = (l^n {I + ikaa)'"- (3.23) 

l^n = V'n (1 + ifcc73)-" (3.24) 

and, then, to consider separately the two columns of these matrices, which we denote as 
follows 

= (Mn>/"n). = {l^n ^ T^n ) ■ (3-25) 

They satisfy the difference equations 

{lTik)tJ.^+i = [l + ikai+Qn]fit (3-26) 

{l±ik)u^^, = [l + ikai+Qn]ut (3.27) 
and have constant asymptotic behaviour 

~ (1 — 1(72) n — > —00 (3.28) 

f„ - (1 - i(72) n +00. (3.29) 



5 



3.2 Existence and analyticity of the Jost solutions fx"^ 

We will show that the Jost solutions fJ-ni^) defined by the summation equations 



E 



(l + ikaiY' 



j = -cc 



(1 - ifc)"-J 



(3.30) 
(3.31) 



are analytic, correspondingly, in the upper half plane and in the lower half plane of the 
spectral parameter k and continuous for fcim > and fcim < 0. 

Let us write the solution of the integral equation (3.30) in the form of a Neumann series 



+ OC 



M+(fc)=E^"W ^nik) 



1=0 



C'j'\k) 

c'/'\k) 



where 



-1 
1 



(3.32) 



(3.33) 



or, in component form, taking into account the identity (3.17) 



1 



2(1 -ifc) 



E 



1 - 



1 + ik 
1-ik 

n-l 



n-j-1 



(ei(''.+-«.)-l)cf2)(fc)| 



(3.34) 



E 



2(1 -ifc) ^ 

j=~oc 



1 



l-i/cj 



l + ^fc \"-J-i' 
1 - ifc / 



(3.35) 



The series ( 3.32| ) is formally a solution of the discrete integral equation ( 3.301 ). One can 
prove by induction on ^ g N that for fcim > 



< 



\i-ikfe\ 



n-l 



_l — — CCi 



J = 1,2. 



(3.36) 



Indeed, for each component j = 1, 2 we have 



< 



1 



2 11 - ik\ 



L-iiei+i-ei) _ ^ 



1— — OG 



l + ik 
1 -ik 



l + ik 
1-ik 



2-1-1 



Cl^^^ik) 
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where we used that the potential 0n is real. Moreover 



1-ik 



-1 




1 + ik 




< 1 + 








1-ik 



< 2 



(3.37) 



since /cim > and rt — ^ — 1 > 0. Therefore we have 



1 I r n 

<^-^T V e-*(«'+-«') - 1 Cf(^'(fc) + Cf'^'\k)] (3.38) 



and we can use the inductive hypothesis to get 

-1 n— 1 



cri'«(fc) 



< 



1 2*^ 



-5(61+1-01) 



l| j ^ |e"i(*'i+i~^'i) - 1 



\li— — oo 



(3.39) 



The use of the summation by parts formula, according to which for any real non-negative 
sequence {bj}^^^ such that the series X^j^oo convergent and for any m G Nq 



fe-i 



fc=— 00 \ j=— 00 



< 



m-\-l 



(3.40) 



completes the proof of (3.36). 
Therefore for fcim > 



< 



P{nf 



J = 1,2 



where 



P(7i) = 2 ^ U(«.+i-e.) _ 1 



(3.41) 



(3.42) 



and we conclude that the Neumann series (3.32) for /i^(A:) is uniformly convergent (with 
respect to n and k in the upper half-plane) and 



provided the potential 9n is such that 



P(oo) = 2 ^ Li{e,+i-e,) _ ^ 



< +00. 



(3.43) 



(3.44) 



j=-oo 

In the same way one can prove that fi~ is analytic for fcim < and continuous for kim < 0. 

3.3 Existence and analyticity of the Jost solutions z/^ 

Let us now consider the summation equations defining the Jost solutions i^^, which can be 
rewritten as 
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1 v-^ (1 — ikai 



^ fi 



1 + (l-ifc)J-"+ 



(1 + ifcV 



(3.45) 
(3.46) 



showing explicitly that their kernel is singular at fc = ±i, i.e. in both the lower and upper 
half planes of fc. This dissymmetry of the summation equations for the Jost solutions 
with respect to /z^ is peculiar of the discrete case, since it disappears in the continuous limit. 
In fact, in order to study the analytic properties of the Jost solutions one needs to define 
the same solutions by introducing alternative summation equations, which can be obtained 
by exploiting the special structure of the potential Qn and the symmetry property of the 
spectral problem (3.1) for n —>■ —n. 
By using the relation 

one can easily verify that 

{l+ik(Tl+Q,,)-^ 



1 



(1 - ik<Ti + cFiQnCri) 



(3.47) 



(3.48) 



and therefore the spectral problem (3T) can be rewritten as 

1 



l + fc2 



If we introduce 



we have 



(1 - ikai + criQ„cri)Xn+l = Xn- 



xn+1 = (1 + k'n 



Cn+1 = (1 - ik(Ji + 0-l(3-„CTi)^„ 

and, consequently, using the same procedure we followed above for </>„ we get 



and coming back to Xn 

Xn^Wn+ J2 (l+fc')""'(l"*fcf^l)'"""'^lQ^-l'^lXf 



(3.49) 



(3.50) 



(3.51) 



(3.52) 



i=n+l 



(3.53) 



where we choose Wn as in ( |3.20| ). 

One can check explicitly that ^„ defined in (3.19) satisfies this summation equation by 
moving in (3.19) the j = n term of the sum from the right side to the left side, by applying 
from the left (1 + ikai + Q„)^"'^(l + ikai) and then by using on the right ( 3.48 ), ( 3.47 ) and 
the spectral equation for ipn. 
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Therefore from (3.53) by considering the transformation (3.24) we get the following 
alternative summation equations for i^^(fc) 

\ / j=,i+i ^ ' 

Following a procedure analogous to that one used for the ^J.^{k) one can prove that the 
v^{k) are also analytic, correspondingly, in the upper half plane and in the lower half plane 
of the spectral parameter k and continuous for /cim > and /cim < provided the potential 



satisfies the condition ( 3.44 ) 



Finally, let us note that, taking into account the explicit expression of the Green's func- 
tions, one can easily obtain the asymptotic behavior of the Jost solutions at large k 

fi+{k) ^ ( ) (fc) ^ ( I ) \k\^oo ki^>0 (3.56) 

fi-ik) ^ f } ) ^^-(fc) ~ ( ) \k\^^ fchn < 0. (3.57) 



3.4 Scattering data 

Let us define the Wronskian of any two vectors v and w as 

W {v, w) = det {v, w) . 



(3.58) 



The vector-valued sequences u„ and Wn are linearly independent if W (vn, Wn) ^ for all n. 

In particular, if v„ and Wn are any two solutions of the scattering problem (B.l), their 
Wronskian satisfies the recursive relation 



W {Vn+l,Wn+l) = (l + fc^) W{Vn,Wn) 

Hence, for any positive integer j 

W{c^t{k),<p-{k)) 



(3.59) 



1 



and 



W{^tik).^nik)) 



-,{k),4>--,{k)) 

{l + k^Yw{^^l_^{k),nl_^{k)) 

;i-f fc2)-^V(v.++^.(fc),^-^^.(fc)) 



^{l + k^Yw{v++Ak),u-^^{k)). 
Taking into account ( [3.28 ) and ( |3.29| ), for k = k^^ in the limit j ^ oo we get 



2\" 



W'(0+(fc),0-(fc)) --2(l + fc2)", M^(V;+(fc),C(fc)) =2(1-Hfc2) 



(3.60) 



(3.61) 



(3.62) 



which shows that the eigenfunctions (j)^ an d 4'~ are linearly independent, as well as ■0^ and 
ip~ . Since the discrete scattering problem (3.1) is a second-order difference equation, there 
are at most two linearly independent solutions for any fixed value of k and consequently we 
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can write 0+ and (j)^ as linear combination of and ■0n or vice- versa. The coefficients of 
this hnear combinations depend on k. Hence the relations (k — /crc) 



^iik)=b^ik)^,^{k) + a^ik)i,^ik) 

which define the scattering data a'^{k), b^{k). 

In terms of the Jost solutions, they can also be written as 



a±(fc) 



p±(fc) 



l±ik 
l^ik 



,,t{k) + i^^{k) 



where we introduced the reflection coefficients 



p^ik) 



a±(fc)' 



Notice that from ( 3.64 ), by using ( 3.2£ ), we can get the asymptotics of /i„ for 



M^(fc)^&±(fc) 



l±ik 
1 =F i/c 



-a±(fc) 



1 



Calculating W ((/)+(fc), 0„ (fc)) using ( 3.63| ) results in 

W (0+(fc), 0;(fc)) = - (a+(fc)a-(fc) - b+{k)b-{k)) W (V'+(fc), CW) 
and taking into account ( [3.62| ) yields 

a+(fc)a-(fc) -6+(/c)6-(fc) = 1 A; e M. 



Moreover, 



or, equivalently. 



W(</.±(fc),V,t(fc)) -T2(l + fc2)"a±(fc) 



a±(/c) = T2W^(M^(fc),^^^(fc)) 



(3.63) 



(3.64) 



(3.65) 



-oo 



(3.66) 



(3.67) 



(3.68) 



(3.69) 



(3.70) 



which proves, from one side, that a~^{k) can be analytically continued in the upper half plane 
(respectively a~ (k) can be analytically continued in the lower half plane) and, from the other 
side, that the zeros of a~^{k) in the upper half plane (respectively of a^(fc)) correspond to 
bound states of the scattering problem. 

We note that for real potentials 6'„, if Xn{k) satisfies the scattering problem (3.1), then 

Xn{k) ^ ta^Xnik*) (3.71) 



satisfies the same equation. Taking into account the boundary conditions ( 3.21 ) and ( 3.22 ), 
we conclude that if the potential 0„ is real the Jost solutions obey the symmetry conditions 



C(fc) = ^<J2 ^nik) = -*^2 (V^+(fc*))' 

Moreover, one can easily show by using ( pS| ) that 

a-{k) = (a+ik*))* 



(3.72) 



(3.73) 
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and by using (3.63) for k ~ k^c that 

6-(fc) = -(6+(A:))*. (3.74) 

The scattering coefficients can also be given as exphcit sum of the Jost solutions. The 
formulae are derived as follows. First, we rewrite for k = /cro the summation equation (3.30) 
as 



^ j — — oo ^ ' 

{1-tk)^ (l + z/fc)J-"+i ^^^^J^"' 



+ 00 



1 + ikai) 



n-j-l 



J=7l 



and, then, using (3.46) we obtain 



t{k)-i^-{k)a+{k)=( M(l-a+(fc))+ £ ii± 



(l + ikai) 



] = -co 
j-n+1 



(1 - ik)"-l ^^^^ 



(3.75) 



Inserting in it (3.64) we have 
5+(fc)^^ + ^^' " 



l~ik 



(l + ifcffl)'' 



(3.77) 



From ( |3.45 ) it follows that the term in square brackets in the left-hand side is (1 ,1)"^ and 
therefore we obtain using (3.17) 



j=-oo 



h+{k) = 



2( 



(fc) 



+ (ei(''^+i-''^^-l)^+'(')(fc) 



(3.78) 



(3.79) 



These (disc rete v ersion of) integral representations, together with the symmetry relations 
(iizl) ^'^'^ ( |3.74|) completely determi ne t he scattering coefficients. 

The discrete scattering problem (3T) can possess discrete eigenvalues (bound states). 
These occur whenever a^(fc^) = for some k'^ in the upper half plane or a~{k^) — for 
some k'[ in the lower half plane. Indeed, for such values of the spectral parameter from 
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( p9| ) it follows that W {ct>+{k+), ^+(k+)) = and IF ((^"(fc^), ^p-{k^)) = and therefore 
the eigenfunctions are linearly dependent, i.e. 



<l^+ik+)=b+^j+{k+) j = l,...,J+ 



(3.80) 
(3.81) 



for some complex constants b^. In terms of the modified Jost solutions ( 3.80|) , ( |3.8l| ) can 
be written as 



1 + ik 



1 — ikp 



Note that when the potential 6'„ is real, (3.73) implies that if fc^ is a zero of a+(fc) in the 
upper half-plane, then kj = (fc^^) is a zero of a^(k) in the lower half-plane and therefore 
J+ = J-. 



4 Inverse Scattering problem 



Let us assume a+(fc) has J+ simple zeros {kj^}^._^ in the upper half plane and a (fc) has 



J simple zeros {kj }''_^ in the lower half plane. Then, by (3.82) it follows 



Res ( ^^^^^■,k+ 



a+{k)''"^ 



da 
~dk 



= /^n {k. ) 



da 
~dk 



Ml- ikt 



= C7 ' 



1 + ifc; 



(4.1) 

(4.2) 



where 



da+ 
Ik 



k=k^ 



da 
~dk 



k=k' 



(4.3) 



Taking into account the analytic properties of the Jo st so lutions and of the scattering co- 
efficients a^(fc), as well as the asymptotic beh avior ( 3.56 ), we can use the Cauchy-Green 
formula to obtain from the "jump" conditions (3.64) 



l + ik+Y ct 



1 
27ri 

1 



k-k: 



C-k + iO \l-iC, 



1 - ifc; 



c. 



ik- 



k~k 



1 
2tH 



3=1 \ ■ 3 / J 

' p-{C) fi-K 



C-k-iO \l + iC 



(4.4) 



(4.5) 
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Equations (4.4) and (4.5) define a Riemann-Hilbert problem on the line which, in principle, 
allows one to solve the inverse problem, i.e. to reconstruct the Jost solutions from the 
scattering data. 

In order to complete the inverse problem, we need to reconstruct the potential from the 
scattering data. Let us rewrite the difference equations (3.27) in matrix form 

[1 



and let us write the asymptotic expansion at large /c of f„ = [v^ , v^) 



AO) 



ik 



(4.6) 



(4.7) 



where, strictly speaking, the expansion of the first column is performed for fcjm > while 
the expansion of the second one is performed for ki^ < 0. Inserting this expansion into (4.6) 



and taking into account that, according to (3.56) and (3.57), 



,(0) 



= 1 — 1(72 yields 



(-1) 



(l + iCT2) 



(4.8) 



5 Time evolution of spectral data 



The time evolution of the Jost solutions and, consequently, of the scattering data is fixed by 



Since the matrix /i„(l + ia^k)" is a general solution of (2.9) there exists a Q{t; k) such 



that 



Xn{t\ k) = ^„(i; fc)(l + lask^nit; k), 



(5.1) 



satisfies both (p. 91 ) and (2.10). 

Inserting (5J) into ( 2.10 ), if 0„ ttt (r G N) for n ±00, by evaluating the limit 
n — 00 we get, recalhng ( 3.28 ), 



(5.2) 



where 77 = (^1)' • 

Then, we evaluate the limit for n +00 by using and we obtain the time evolution 

of spectral data (fc = fcRo) 



btik) = TTr^^b^{k), at{k) = 
k =F 107771 



(5.3) 



I.e. 



6±(fc,i) = &±(fc,0)e=F2*T*, a±(fc,t) =a±(fc,0). (5.4) 
For the discrete spectral data one gets that the locations of poles at k = k^ are fixed while 



(5.5) 



Notice that the switching on of the time evolution introduces a singularity at fc = 0. 
One can prove that this does not spoil the good properties of the spectral transform. In the 
following section where the matrix spectral problem is reduced to a scalar problem this will 
be shown explicitly. 
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6 Lax triplet 



6.1 Semi-discrete case 

We can write the spectral problem (^^) in component-form and decouple the two equations 
to get 



An+2 ~ ynXn+l (,-1- ^ «- )Xn 
1 



.(2) 



ik 



An+l 



e 2^ 



+ 1-9^)^(1) 
An 



(6.1) 
(6.2) 



with 



(6.3) 



i.e. we recover the discrete Schrodinger equation whose scattering theory was given in 
In this case the time evolution is fixed by the pair of spectral problems 



(6.4) 
(6.5) 



The compatibility requirement for the three spectral problems (|6.1| ), ( |6.4[ ) and (6.5) furnishes 
after one integration the semi-discrete sine-Gordon equation ( |2.8| ). 
Note that, if we introduce 



the equation defining can be rewritten as 

gn = Otn+l 



(6.6) 



(6.7) 



that can be considered as a discrete Riccati equation, which can be solved for in terms 
of a continuous fraction of g„ (see pO[). 

Note also that for 0„ vanishing at large n (modvr) we have from (6.3) that — 2 — + 



and the spectral theory developed in |9|] is applicable. However, the condition required to 
be satisfied by (7„ — 2 in ||] is stricter than (3.44). 

We conclude that the Orfanidis matrix spectral problem can be reduced to a scalar 
Schrodinger spectral problem if the potential On (mod tt) sufficiently rapidly. 

6.2 Time evolution of spectral data 

The evolution of the spectral data can be determined by considering the first and the third 



Lax operators, i.e. the spectral problem ( |6.lD and (6.4). 

In order to meet the notations used in M we introduce the eigenfunction as follows 



The principal spectral problem ( |6.lD reads 

(1 + ik)^ri+2 + (1 - ik)£.n = gniri+l 



(6.8) 



(6.9) 
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while the auxihary spectral problem (6.4) becomes 



(6.10) 



We need the Jost solution /i^(fc) introduced in which is analytic in the upper half- plane 
and is defined via the following discrete integral equation 



1 



l + ik 
1-ik 



(5,_i-2)/i+(fc) 



j=n+l 

and which has the following asymptotic behaviour for n +oo 



lim ^l+ik)^l, ki^>0 



and for 



n — (X) 



fit{k)-a+{k) + 



1-ik 
l + ik 



b+{k), kin, = 



(6.11) 



(6.12) 



(6.13) 



where a~^{k) is the inverse of the transmission coefficient and 6+ (A:) is the reflection coeffi- 
cient. 

Then we look for a solution of ( 6.10| ) of the form 

Ut;k)^nit;k)^l+{t■,k) (6.i4) 

If we substitute ( 6.14 ) into ( |6.10 ), we get 

e [nt^l~'^^+ + = 7e-^(''"+^+''"V+ - (l + «fc)7e-^'>++i. (6.15) 

If 9n — *■ rTT (r e N) for n — > ±oo, taking into account the asymptotic behaviour of /j,^ as 
n +0O we get, first, 



k + iOrj^t 



V 



-ly 



(6.16) 



and, then, by considering the limit for n —^ — oo {ki^ ~ 0) the evolution equation for the 
spectral data 



a+ik) = 



btik) 



k — iOrjjt 
which can be trivially integrated to 

a{t;k) = a{0;k) 

h{t; k) = 6(0; k) exp 



b+{k), 



2ir]j 



(6.17) 
(6.18) 

(6.19) 
(6.20) 



Let us now study in details the behaviour of the Jost solutions at A: = when the time 
is switched on. For the sake of simplicity, let us take r]j = 1, i.e. let us scale the time. 
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In 1^ the Ricmann-Hilbert boundary value problem defining the Jost solutions was given 

as 



a+(fc) 27ri \l + is/ s — fc — zO 

For the sake of simplicity, we omit a possible contribution from the discrete part of the 
spectrum. It can be added without difficulty along the same lines followed for the Orfanidis 
spectral problem. 

For convenience we introduce 



Deriving ( 6.21 ) with respect to time we have 



--^l^Vi^U^^tfMH^,^ (6.24) 
l + isj s-k + iO ^ ' 



and therefore 



k + iOt 2m 7_oo V 1 + s-k + iQ \s ~ iOt k + iOt 

1 fl-isY <fr,{-s)p+{s) 



Stti J-oc \^ + is J s — k + iO 

Since 

1 1 



s-iOt k + iOt [k + iQt) [s - iOt) 



(6.25) 



(6.26) 



we have finally 



fc + iOi 27ri 7_oo + s — fc + iO 
where 

5„(0 ^ 1 + ^ f 1^)" A^^izfV^if) (Si for . > 3 

27rii_^ Vl + «V s-iOt \ p-{0) fori<0 

In order to study the singularity of <&„(fc) at fc = one can try to solve the integral equation 
by iteration and see at each step how the singularity 1/ (k — iOt) is transformed. 

Let us first consider the case t > 0. Iterating once yields an integral equation that 
contains the distribution 



11 1/1 1 



iOs- k + iO k~iO\s~k + iO s ~ iO 



(6.29) 
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Therefore in this case the iteration renormahzes the coefficient of the singularity 1/ (fc — iOt) 
but does not change the nature of the singularity at fc = 0. 

Let us now consider the case t < 0. In this case in iterating we get in the integral 



1 



1 



s + iOs- k + iO 



(6.30) 



which is a distribution continuous at A; = 0. Therefore in this case we can conclude that the 
singularity of $„(A;) can be singled out and it is just 



k - m 



(6.31) 



Recalling the definition of $„ we conclude that, for t < 0, dt^~{k), thanks to the continuity 
of fJ-nik) at fc = 0, is less singular than 1/fc. 

In conclusion, fc9t/i^(fc) and, therefore, k'^dt£,~{k) appearing in the auxiliary spectral 
problem, are continuous at fc = 0. Analogously of course for fc9t/i+(fc) and k'^dt(,^{k). They 
are, however, discontinuous sd t — 0, i.e. g„(t) and, consequently, On{t) evolve according 
to different laws for t ^ 0, which is not surprising since the sine-Gordon equation is not an 
evolution equation. 



6.3 Doubly discrete case 



Analogously in the system of the two spectral problems (2.1) and (2.2) the two components 
Xn)m and Xriln of Xn,m Can be decoupled. Of course Xn)m satisfy the same discrete version 
of the Schrodinger spectral problem as Xn \ i.e. 



(1) 



^ gn,mXn+l,m - (1 + ^ )X 

where 

„ _ „5(ere+l.m-en,n,) _|_ „ - 5 (6^ + 2 , m " f „+ 1 , „ ) 

yn.ni — ^ c , 

while the time evolution is fixed by the couple of spectral problems 



3.32) 



(6.33) 



An,m+ 



l = {k 



76 



^(1) ^ y(l) 



76 



.+e„.„+i) 



i+e„ 



^(1) 



76 



-5(e„,„+i+e 



An,m+1 



y(l) 



An+l,m 



(6.34) 



(6.35) 



The compatibility among the three spectral problems ( 6.32 ), ( 6.34 ) and ( 3.35 ) furnishes, 
after an integration, the Hirota-Bianchi doubly discrete sine-Gordon equation. 



6.4 Time evolution of spectral data 

Following a procedure analogous to that followed in the semi-discrete case we look for a 
solution of ( |6.34| ) of the form 



^nMk) = nrn{k)fi+{k){l+iky 



(6.36) 
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where /i+„(/c) is the Jost sohition introduced in Taking into account the asymptotic 
behaviours ( 6.12| ) and (3.13), if 9„.rn rn {r € N) for n ±00, we obtain, respectively, for 



n +00 the time evolution of ft 



a^+i(fc)f}7„i(fc) = l-^ (6.37) 



and for n —>■ —00 the time evolution of the spectral data (ki^ = 0) 



am+i{k) ^ a„i{k) (6.38) 
(777 - ik) ^ 
(777 + ik) 



Wi(fc) = -P^-^braik) (6.39) 



where 77 is defined as above. 
If we introduce 



a^n = CL{t + Tin) , bjn = b{t + T'm), t ^ nr, (6.40) 
change the constant 7 as follows 

7 ^ T7 (6.41) 

and take the limit t and n 00 at t fixed we recover, as expected, the time evolution 
of the spectral data for the semi-discrete sine-Gordon equation. 
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